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ABSTRACT 

Aims. We investigate the distribution and properties of Faraday rotating and synchrotron emitting regions in the Galactic ISM in the 
direction of the Galactic anti-centre. 

Methods. We apply Faraday tomography to a radio polarization dataset that we obtained with the WSRT. We developed a new 
method to calculate a linear fit to periodic data, which we use to determine rotation measures from our polarization angle data. From 
simulations of a Faraday screen -I- noise we could determine how compatible the data are with Faraday screens. 
Results. An unexpectedly large fraction of 14% of the lines-of-sight in our dataset show an unresolved main component in the Faraday 
depth spectrum. For lines-of-sight with a single unresolved component we demonstrate that a Faraday screen in front of a synchrotron 
emitting region that contains a turbulent magnetic field component can explain the data. 

Key words, magnetic fields - radio continuum:ISM - ISM: magnetic fields, structure - techniques: polarimetric, miscellaneous 
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1. Introduction 

Faraday tomography is a very powerful tool for studying the rel- 
ative line-of-sight distribution of regions with synchrotron emis- 
sion and Faraday rotation. If the thermal plasma does not emit 
its own synchrotron radiation it simply acts as a Faraday rotat- 
ing screen for any linearly polarized radiation that illuminates 
it from the back. However, if the thermal plasma co-exists with 
relativistic particles which emit (linearly polarized) synchrotron 
radiation, the situation can be much more complex. This is be- 
cause the contributions from various parts of the line-of-sight 
through the plasma combine vectorially, which may even lead to 
complete cancellation (at certain wavelengths) of the polarized 
signal emitted in the plasma. 

For a simple Faraday screen, the amount of rotation of 
the plane of linear polarization gives direct information on the 

observer 

Faraday depth 'R [rad/m^] of the screen: •R = 0.81 / «,B-d/, 

source 

the line-of-sight integral of the product of the electron density 
[cm^^*] and the magnetic field component [/iG] along the line- 
of-sight [pc]. In that case, fi and the rotation measure RM - 
d^ldA^ (where O is the polarization angle of the radiation) 
are identical. In the more general case, with a mixture of syn- 
chrotron emitting and Faraday rotating layers along the line-of- 
sight, the relation between RM and R becomes more compli- 
cated. SokolofF et al. (1998) have calculated the result for var- 
ious geometries of synchrotron emitting and Faraday rotating 
layers. When the synchrotron emissivity and Faraday rotation 
are both constant per unit line-of-sight (a so-called 'Burn slab'. 
Burn 1966 ), RM - 0.5'Rmax, where 'Rmax is the Faraday depth of 
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the far side of the slab. Sokoloff et al. showed that this is true for 
any distribution that is mirror-symmetric along the line of sight. 

We present a first discussion of the use of WSRT observa- 
tions at A ~ 0.86 m, of the linearly polarized component of the 
diffuse Galactic synchrotron emission towards the Galactic anti- 
centre. By using Faraday tomography we investigate the proper- 
ties of the magneto-ionic ISM in this direction, and in particular 
the abundance of Faraday screens. 



2. The data 

The present dataset was obtained with the WSRT, a 14-element 
E-W interferometer of which 4 elements are moveable to im- 
prove (u,v) coverage. We used the mosaicking technique to map 
an area of about 7x7 degrees^ with 49 pointings. The distance 
between pointings was chosen such that instrumental polariza- 
tion is suppressed to less than 1% (Wieringa et al. |1993l l. In our 
analysis we leave out the edge of the mosaic where instrumen- 
tal polarization effects cannot be suppressed by mosaicking. We 
also exclude lines-of-sight for which instrumental polarization 
of off-axis sources is an important factor The central coordi- 
nates of this field (in the constellation Gemini) are a - 7'' 18™ 
and 6 = 36°24' (J2000.0), which is Z ^ 181° and ^ 20° 
in Galactic coordinates. The observations cover the frequency 
range between 315 and 385 MHz, with 213 usable independent 
spectral channels of about 0.4 MHz each. The field was observed 
for 6 nights (@ 12 hrs each) in December 2002 and January 
2003. This yielded visibilities at baselines from 36 to 2760 me- 
ters, with an increment of 12 meters. The results we discuss here 
were obtained with a Gaussian taper that decreases the resolution 
to 2.2' X 3.7' (FWHM at 350 MHz). 
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Fig. 1. <1>(/}^) distribution (top panel) and spectrum (bottom 
panel) for a line-of-sight with a main peak that is similar to a 
Faraday screen (low A). The of the RM fit is 1.02 (deter- 
mined by using the RM with the maximum P(RM) - see the text). 



The data were reduced using the NEWSTAR data reduction 
package. Dipole gains and phases and leakage corrections were 
determined using the unpolarized calibrators 3C48, 3C147 and 
3C295. Due to an a-priori unknown phase offset between the 
horizontal and vertical dipoles, signal can leak from Stokes U 
into Stokes V. We corrected for this by rotating the polariza- 
tion vector in the Stokes (U,V) plane back to the U axis, assum- 
ing that there is no signal in V. The polarized calibrator sources 
3C345 and DA240 defined the sense of derotation (i.e. to the 
positive or negative U-axis). Special care was taken to avoid au- 
tomatic flagging of real signal on the shortest baselines. From 
2 lines-of-sight with a strong polarized signal we estimate that 
the amounts of ionospheric Faraday rotation in the 6 nights are 
identical to within ~ 10° so we did not correct for that. 

Structure on large angular scales in Stokes Q and U will have 
been filtered out because we have no information on baselines 
below 36 m. We discuss the influence of this on our results in 
Sect. [3:21 



3. Analysis 

3.1. Methods 

The wide coverage in space (from 0.6m^ ^ A^ ^ 0.9m^), 
combined with the relatively small channelwidth, allows us to do 
Faraday tomography, also known as Rotation Measure Synthesis 
(see e.g. Brentjens & De Bruyn I2005I I. which probes the dis- 
tribution of Faraday rotating and synchrotron emitting regions 
along the line-of-sight. In Faraday tomography the polariza- 
tion vectors of the individual channels are 'coherently added' 
by derotating the vectors using an assumed Faraday depth R: 
PCR) = J P(A^) e^''^'-'^%-^'^^'dA^, where P(A) and (D(/l) are the 
polarized intensity and polarization angle of an individual chan- 
nel at wavelength A. PCR) - IPCR)] is the intensity of the po- 
larized emission at Faraday depth R., and a Faraday depth spec- 
trum (or R spectrum) can be constructed by calculating P{R) and 
arg{P{R)) for many values of R. We calculated R spectra be- 
tween -72 rad/m^ and +72 rad/m^ in steps of 1 rad/m^. Previous 
surveys indicate that the R of the diffuse emission in this volume 
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Fig. 2. Identical to Fig.[T] However, here the;^^^^ of the RM fit is 
10.58, which is due to a complicated R spectrum. 



of the Galaxy lie well within this range (see e.g. Spoelstra |1984| l. 
The Point Spread Function (PSF) of the f ('R)-determination, 
which is the Fourier Transform of the /1^-sampling of the data, 
has for our data a FWHP of about 12 rad/m^. From simulations 
of the noise in the individual channels we estimate that 99% of 
the noise realisations will lie below the 0.7 K level in our "R spec- 
ti-a. 

We quantified the behaviour of the H spectrum in two ways. 

First, we derived a measure of the symmetry of the R spec- 
trum as a whole, by calculating the reduced Ck'red-' ^ linear 
fit to the <l)(/l^) data. A Faraday screen or a symmetric distribu- 
tion of Faraday rotating and synchrotron emitting regions along 
the line-of-sight (for example a Burn slab) will show a linear 
<1>(/1^) dependence. The;t'^^j of a linear fit to the C)(/}^) distribu- 
tion can therefore separate symmetric and more complex distri- 
butions along the line-of-sight. 

Fitting periodic data like polarization angles requires that 
the nl80° periodicity of 0(/l) is properly taken into account. In 
Schnitzeler et al. (2007) we introduced a method that finds the 
lowest-^j^gj straight-line fit to polarization angle data by going 
through all the possible 180° wraps of polarization angles that 
are allowed by the data. However, the number of configurations 
of wraps of the individual datapoints increases strongly with the 
number of datapoints. For the present dataset, with of order 200 
<l>(/l)-values for each line of sight, application of that method is 
thus not practical. 

However, by writing <l)(/l^) as the complex number e^'*^'''', 
the RM spectrum can be calculated: P{RM) - \P{RM)\ - 
I J ^2mA^)^-2iRMA^^;^2y j^iQ powcT P{RM) shows which RM 'fre- 
quencies' create the observed ^(A^) dependence. If O depends 
linearly on A^, the RM with the maximum P(RM) will be the 
best fitting slope for the data. If the RM spectrum is more com- 
plicated, the RM with the maximum P(RM) will not give the 
best linear fit to the ^{A?), and the x^f.^ of this fit will be higher 
than that of the best linear fit to the (i>{A^). Even the best linear 
fit to the <l)(/l^) will then have a large 

The second criterion quantifies the deviation of the main 
peak in the R spectrum from a Faraday screen, in the presence 
of noise. We define A as the root-mean-square vertical separa- 
tion in the R spectrum between a peak in the R spectrum and the 
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Fig. 3. Distribution of independent lines-of-sight with a main 
peak that satisfies the A criterion (green pixels), or withx^^,^ < 2 
and max(P('R)) > 2 K (yellow pixels), or that satisfy all three 
criteria (red pixels). Lines-of-sight that did not pass these crite- 
ria are shown as black pixels. The pixel size is 2.2'x2.2'csc(5). 
The outer edges of the mosaic show an increase in x^^^, which 
explains the absence of yellow pixels. 



PSF that is scaled to the same height as the peak. A is calculated 
over the PSF out to the point where the PSF goes through its 
first minimum, at a distance of 10 rad/m^ from the centre of the 
peak. A Faraday screen will show up as an unresolved peak in 
the spectrum, which means that its A will be lower than that 
of a peak in the K spectrum that is too broad to be fitted by a 
PSF. By comparing the A we find for the main peak in each line- 
of-sight to the distribution of A of a Faraday screen + noise, we 
can quantify how improbable it is that the main peak is due to a 
Faraday screen. 

In Figs.[T]and|2]we show the <t>(A^) distributions and spec- 
tra for two lines-of-sight that have a main peak that is similar 
to a Faraday screen (i.e. low A). One of the lines-of-sight has a 
simple K spectrum and low;^^^^^, the other has a complicated K 
spectrum and high;^^^^^. Fig. |2] shows that low A are also found 
for lines-of-sight with higher In this case the high;^'^^^ is 
due to the presence of multiple components in the spectrum. 

3.2. Results 

From a total of 22.800 independent lines-of-sight in our dataset, 
1757 lines-of-sight have a x^g^i < 2 and a maximum in the 'R 
spectrum > 2 K, about 2.5 times the level below which 99% 
of the noise PCR) lie. In Fig. [3] we plot the sky distribution 
of these 1757 lines-of-sight as yellow pixels. We indicate the 
lines-of-sight that furthermore have a main peak whose A is 
less than 99% of the A that we found in our simulations of 
a Faraday screen + noise as red pixels. The pixels in Figs. 
[3] to |5] indicate independent lines-of-sight. These figures cover 
Galactic longitudes from (Z, Z7)=(178°,26°) in the top left cor- 
ner to (l, Z7)=(185°,15°) in the bottom right corner. The Galactic 
plane has a position angle (north through east) of about 21° rel- 
ative to the vertical axis in these plots. 
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Fig. 4. PiK) of the highest peak in the R spectrum for lines-of- 
sight where the secondary peak in the R spectrum has at most 
half the strength of the main peak. Lines-of-sight that do not 
pass the selection criteria are shown in black, together with 455 
lines-of-sight that have PCR) < 2 K. 



In Figs. ID and |5] we plot the distribution on the sky of PCR) 
and fi of the main peak in the fi spectrum. Ionosphere mod- 
els indicate that it contributes only 0.6 rad/m^ to the "R shown 
in Fig.|5](Johnston-Hollitt, priv. comm). For these maps we se- 
lected lines-of-sight where the 'R spectrum is dominated by a 
main peak, by requiring that this peak should be at least twice 
as high as the second brightest peak. Lines-of-sight selected in 
this way show a tight coiTelation between and RM, with a Icr 
scatter of about 0.2 rad/m^. 

Structure on large angular scales in Stokes Q and U (and 
therefore in polarized intensity P and polarization angle O) will 
not be picked up by an interferometer. However, if large-scale 
structure were produced at the same K as the emission visible in 
Fig. m the large gradients per pointing in Fig. |5] create struc- 
ture in O on small enough angular scales that can be picked up 
by the WSRT. Large-scale structure would then become 'visible' 
by passing through the Faraday rotating foreground with these 
gradients, and Figs. |4] and |5] therefore form more or less 'com- 
plete' maps in that they contain most of the emission at those 
Faraday depths. 

3.3. Faraday screens 

The main peak in 3086 lines-of-sight (14% of the total number of 
independent lines-of-sight) has a A that is less than 99% of the A 
we found in our simulations of a Faraday screen + noise, making 
these main peaks indistinguishable from Faraday screens. These 
lines-of-sight are indicated as green pixels in Fig.|3] These peaks 
also closely follow the H - RM relation expected for Faraday 
screens, with a Icr scatter of about 0.2 rad/m^. The actual frac- 
tion of lines-of-sight with Faraday screen components in their H 
spectrum can be higher than 14% since we only tested the main 
peak of the "R spectrum with our A criterion. 

Radio frequency interference and/or an underestimation of 
the noise levels will increase the values of A expected for a 
Faraday screen -i- noise. The lines-of-sight that already pass this 
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Fig. 5. fi corresponding to the PCR) plotted in Fig.|4] The lines- 
of-sight plotted in these figures are selected by using the same 
criteria. Lines-of-sight indicated in white have K > 30 rad/m^. 
Lines-of-sight indicated in black either did not pass the selec- 
tion criteria or have < -6 rad/m^ (652 and 9734 lines-of-sight 
resp.). 



criterion with our conservative choice for calculating A will then 
still pass the A criterion; there could however be more lines-of- 
sight that would qualify as Faraday screens. 



4. Discussion 

We first estimate the implied value of the large-scale magnetic 
field, for a representative value of 7? - +10 rad/m^ (Fig.|4|. If 
the Faraday rotating region has a depth L and consists of clumps 
of electrons with a constant density that occupy a fraction 
/ of L, with n,=0.08 cm"^ and /=0.4 (Revnolds |l991l l. L fol- 
lows from observations of the emission measure, EM. The EM 
from the Wisconsin Ha mapper (WHAM; Haffner et al. I2003I I 
are on average 2 Rayleigh in this region, or 4 cm^^pc for a 10"^ K 
gas, which indicates a line-of-sight of 1 .5 kpc length. For these 
lines-of-sight extinction of the Ha line plays only a minor role. 
The DM-50 cm^^'pc implied by this model agrees well with the 
DM=54 cm"^pc calculated from the NE2001 model by Cordes 
& Lazio (I2003I I for a line-of-sight in the same direction and of 
the same length. 

For this model the integral for is also easy to solve. By 
using the same and L, the large-scale magnetic field com- 
ponent along the line-of-sight, Breg,||, must be about 0.2 yuG to 
explain the average "R-IQ rad/m^ present in our data. If part of 
the Ha emission is coming from beyond the Faraday rotating re- 
gion, we would overestimate the depth L of the Faraday rotating 
layer, and the Breg.n we derive would be too low. A reasonable up- 
per limit for Bieg.n can be found by adopting a pitch angle of 8° 
and a total strength of 4 fiG for the large-scale field (both values 
from Beckl2007l), which gives Bieg,||,max ~ 0.6 fiG. Note that the 
change of sign of 'R in Fig.|5]indicates that the observed features 
are also partly due to variations in the magnetic field geometry 
and not only to variations in the electron density. 

From the Haslam et al. ( 119821 ) 408 MHz data we estimate 
that the Galactic foreground in this direction has a brightness 



temperature of 36 K at 350 MHz, assuming a -2.7 spectral in- 
dex for the synchrotron brightness temperature. For a completely 
uniform magnetic field and without internal Faraday rotation the 
peak in the R spectrum would have a polarized brightness tem- 
perature, Tb.poi, of 70% of this, or 25 K, but the Tb.poi of the peak 
in the R spectrum is on average only 5 K. A similar discrepancy 
between the expected and measured Tb.poi was noted for a field 
observed with the WSRT at (l,b)^{llQ°,ir) (De Bruyn et al. 
f2006] l. 

If the R distribution has a finite width, the 25 K is diluted 
over the entire width, and the peak of the distribution will have 
lower PCR). Any further reduction of PCR) requires a turbulent 
magnetic field in the emitting region. For lines of sight with a 
single, essentially unresolved peak in the fi spectrum (i.e. with 
a low value of A) it is not trivial to separate the two effects. We 
simulated a layer with both uniform and Gaussian distributions 
of coexisting Faraday rotation and synchrotron emission. From 
these simulations we conclude that it is not possible to reduce 
the Tb.poi of the maximum in the PCR) spectrum below 20 K, 
and be consistent with the observed low value of A. Therefore 
a turbulent magnetic field component is required to explain the 
observed Tb.poi. 

Another possible picture for these lines of sight with small A 
could be a combination of an !R-extended structure and an emis- 
sion region with a turbulent magnetic field. However, since the 
•R-extended structure has to be narrow to pass our A criterion, 
the depth of this structure will not be enough to explain the 10 
rad/m^ average Faraday depth in Fig.|4] A Faraday screen is then 
required to produce the bulk of this Faraday depth. 

An extreme example of this would be an emission region 
with a turbulent field, observed through a Faraday screen, which 
would be the simplest and most straightforward explanation. In 
that case we can also estimate the relative strengths of the turbu- 
lent magnetic field and the magnetic field component perpendic- 
ular to the line-of-sight, Z?turb and Breg.j. resp. If the synchrotron 
emitting region contains both Biee,± and Z?turb, the expected po- 
larization fraction is reduced by a factor of B^^^ ±/iB^eg ± + ^furb-* 
(Burn'1966V To reduce the expected 25 K polarized brightness 
temperature to 5 K then requires Bturb /^reg.j. ~ 2. This ratio is an 
upper limit if there are other components in the R spectrum that 
contribute to the predicted Tb^poi = 25 K. 

The region around (a, 5)=109°,35°shows some conspicuous 
features. It is bright in polarized intensity, shows a low x^^^, and 
also contains a number of lines-of-sight that could be Faraday 
screens. The band between 36° < ^ < 38° contains many lines- 
of-sight with high x^g^ and broad main peaks and/or significant 
secondary peaks in the R spectrum. We will discuss these regions 
in a future article. 

5. Conclusions 

We used high spectral resolution radio polarization data to study 
the Galactic ISM, and we demonstrated a number of tools that 
can be used for this purpose. In the present dataset we identi- 
fied a significant number of lines-of-sight that show an unre- 
solved peak in their Faraday depth spectrum, similar to a Faraday 
screen. We also studied the spatial behaviour of the principal 
component in the Faraday depth spectrum, and for the lines-of- 
sight that only have one unresolved component in their H spec- 
trum, we showed that a model of a Faraday rotating region in 
front of a synchrotron emitting region that contains a turbulent 
magnetic field component can explain the data. 
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